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EFFECTS OF TEMPERATURE -PROFILE VARIATION 
ON REFRACTION OF SOUND BY JET FLOW 

By Jay C. Hardin 
Langley Research Center 

SUMMARY 

In this analysis, the techniques of geometrical acoustics are employed to study the 
effect of variation of the jet-temperature profile upon propagation of sound in the jet. It 
is shown that the temperature profile of the jet can be shaped to cause the jet flow to act 
as a channel for the sound. Further, a simple analytical expression which predicts the 
occurrence of such phenomena in a nonspreading jet is obtained. This expression also 
seems to be approximately valid for a spreading jet. 

INTRODUCTION 

The phenomenon of refraction of sound by temperature and velocity gradients in the 
medium through which the sound propagates has long been recognized. (See ref. 1.) In 
certain instances, these gradients form ducts which can trap and channel the sound over 
large distances. (See ref. 2.) Because of considerable interest in the noise of aerody- 
namic jets, this result suggests that modification of the velocity and/or temperature pro- 
files of a jet might favorably alter the directivity pattern of the sound without increasing 
its intensity. 

A recent experimental study (ref. 3) has shown that the propagation of sound in a 
jet flow is, in fact, highly dependent upon the temperature and velocity of the jet. Fur- 
ther, while the intensity of the jet noise has been found to be a function of both the magni- 
tude and shear of the mean velocity (ref. 4), it has also been found to be practically inde- 
pendent of temperature (refs. 5 and 6). Thus, it may be possible to achieve beneficial 
alterations in the jet -temperature profile without increasing the production of noise by 
the jet. 

An experimental study in which such refraction was observed is reported by 
Grande (ref. 7). In this study a very cold nitrogen jet was used to study the propagation 
of sound from a source placed in the jet as well as from the jet itself. The results indi- 
cate that the sound from both the source and the jet is refracted inward and obtains a 
maximum intensity along the jet axis instead of the heart-shaped directivity pattern 
normally obtained in jet flows. 


In this analysis, the propagation of sound in a jet flow is examined. The wave 
equation for a finite, spreading jet is complex (ref. 8) and has not yet yielded usable 
results. For this reason, the techniques of geometric acoustics have been employed, 
although they are not strictly applicable. This inapplicability is due to the fact that the 
use of such techniques requires the gradient of the speed of sound to be small with 
respect to the ratio of the speed of sound to the wavelength of the sound (ref. 9). This 
condition occurs only for the high-frequency components of the jet noise. However, 
increasing the frequency of the sound introduced into a jet causes an outward rotation of 
the peak of the directivity pattern (ref. 3). Thus, if it can be shown that the high- 
frequency components can be channeled, there is some reason to believe that the lower 
frequency components may also be channeled. 

SYMBOLS 


A 

constant in equation (27) 


a 

speed of sound 


a a 

ambient speed of sound 
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b i 

constants, where i = 1,2,3 

(eq. (30)) 

c i 

constants, where i= 1,2,3 

(eq. (31)) 

d 

jet diameter 


F 

function defined by equation (12) 

f(y) 

function defined by equation (19) 

T,T, k 

unit coordinate vectors 


Z,m,n 

direction cosines 


M a 

ambient Mach number 
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R 


universal gas constant 


r radial coordinate 

T a ambient temperature 

T H high -temperature -core profile 

Tj jet temperature 

low -temperature -core profile 
T g space average temperature deviation 

t time 

2a n t 
t' = — — 
d 

u component of velocity along X-axis 

u velocity vector 



v disturbance velocity vector 

X,Y,Z jet axes 

x,y,z Cartesian coordinate system 

x 0 ,y 0 initial position coordinates 

x\y',z' nondimensional Cartesian coordinate system 

y c singular value y at which ray becomes horizontal 
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y* constant in equation (28) 

a constant which determines shape of profile (eq. (24)) 

/3 constant which determines rate at which jet spreads (eq. (24)) 

y constant (see eq. (27)) 

constant (see eq. (27)) 

62 constant (see eq. (28)) 

6 angle defined in figure 2 

rj unit vector normal to the wave front 

X c constant (see eq. (18)) 

p adiabatic constant 

<7 spread variable 

<p spread angle 

THEORETICAL ANALYSIS 

Equations of Geometric Acoustics for a Spreading Jet 

A cylindrical jet with the geometry shown in figure 1 is considered. The jet is 
assumed to have a velocity profile 

u = u(x,r) T (1) 

where 

r = (y 2 + z 2 ) 1//2 

which neglects the small radial component of the velocity field and any turbulence which 
may arise in the jet. Likewise, the jet temperature is given by 
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Tj = Tj(x,r) 


( 2 ) 


in degrees absolute. Thus, the speed of sound at any point in the jet is 


a(x,r) = jp,RTj(x,rjj ] 


where R is the universal gas constant and fJ. is the adiabatic constant. 

An acoustic disturbance which arises in the jet propagates outward from its source 
with a ray velocity, obtained from reference 10, of 

v = arj + u (4) 

where rj = zT + mX + nk is a unit vector normal to the wave front. The change in this 
vector along a ray, given in reference 10, is 


^ = va + (vu) • r? - rj(rj • Va + Jj • (vil) • rjj 


For the velocity given in equation (1) and the sound speed given in equation (3), equa- 
tion (4) yields 

dx , 

~ = u + la. 


^ = ma 
dt 


and equation (5) yields 




dm + i §u + 9a = mF 
dt 9y 9y 


dn + i 9u + 9a = nF 

dt 8z 8z 


5 



where 


< 12 > 

Now, a ray which starts at the origin of the coordinate system is considered. By 
treating the ray and the X-axis as vectors and by noting that any two intersecting vectors 
define a plane, the plane defined by these two vectors may be assumed, with complete 
generality, to be the XY-plane by a simple rotation of the coordinate system. Then 

Z(0) = cos 9 m(0) = sin 9 n(0) = 0 

where the angle 9 is defined in figure 2. Thus, equations (8) and (11) may immediately 
be integrated to yield 


z(t) = 0 n(t) = 0 

Therefore, a ray which starts in a plane will remain in that plane, and the three- 
dimensional problem has been reduced to two dimensions. The governing equations 
become equation (6), 


equation (7), 


dx 

dt 


u + Za 


$L = 

dt 


ma 


and 


dZ 

dt 


dm 

dt 



9u l m 9a ,2 9u 
9x 9x ” 9y 


^+m^-Z 2 9u 
9x 9x 9y 



with the initial conditions 


(13) 

(14) 


x (0) = y(0) = 0 


Z (0) = cos 9 


m(0) = sin 9 
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This set of coupled, nonlinear, partial differential equations is complex but can readily 
be solved numerically. Examples of such solutions are given in a later section. 

Channeling in an Axially Uniform Jet 

Some insight into the ray solutions is gained in an analytical sense by neglecting the 
axial or X-dependence of the velocity and temperature profiles. Equations (6) and (7) 
remain unchanged 


dx 

dt 


= u + 


l a 


& = 
dt 


ma 


and equations (13) and (14) become 



dm _ , 2/7 du da\ 

dt ~ L ^dy + dyj 


When m ^ 0, equations (7) and (15) combine to yield 

<M 1 da ^ _ 1 du -,2 
dy a dy a dy 


(15) 

(16) 


(17) 


Equation (17) is known as a Bernoulli differential equation (see ref. 11) and admits the 
solution 


K y) = 


^c 



(18) 


where, by the initial conditions, 


= _a(0)_ 


cos 6 


+ u(0) 


If for all y, -1 < y) < 1, then since + m^ =1, m / 0 for any y and the solu- 
tion (18) is valid everywhere. In addition, since m = sin d, the condition that m ^ 0 
implies that the ray never becomes horizontal and cannot be trapped in the jet. This 
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implication can also be seen from equation (7). Since m / 0, dy/dt never vanishes 
and y(t) must be monotonically increasing or decreasing. Thus, for a ray to be 
channeled, l( y), given by equation (18), must be equal to ±1 at some y. 

Suppose that at y = y c , the ray becomes horizontal. Then, since the velocity and 
temperature profiles are even functions of y, the ray will also become horizontal at 
y = -y c . Then from equations (6) and (7), the equation for the rays becomes 

( 1 ! 

W (u + ia) 2 


If a(y) and u(y) are given, then f(y) is completely determined. The class of equa- 
tions exemplified by equation (19) has been studied extensively and admits solutions in 
terms of elementary, elliptic, or hyperelliptic functions, depending upon the power of y 
in f(y). (See ref. 11.) However, for general f(y), certain statements may be made 
about the nature of the integral curves. (See ref. 12.) One statement which is relevant 
to this discussion follows: 

If an integral curve starts in the fundamental strip -y c < y < y c , it cannot pass out 
of that strip. Thus, channeling will occur for this curve. In fact, the initial condition 
that the ray start at the origin of coordinates may be relaxed somewhat to yield the gen- 
eral principle: 

Let there be given a velocity u(y) and sound speed a(y) which are even functions 
of y. Then if a ray starts at any point ^x 0 ,y Q j with initial normal vector in the XY -plane, 

a necessary and sufficient condition for the ray to be channeled is that there exists a 
y c > 0 such that y c > y Q > -y c and 


l 




±1 


where X c = c ^ s ^ g- + u (y Q ) and cos 0 is *h e x-component of the initial normal vector. 

It might be noted that if the initial normal vector does not lie in the XY-plane, the 
two-dimensional character of this analysis does not apply, and the rays would be helical 
in nature. 
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NUMERICAL ANALYSIS 


2x 2v ^ u * ^^3 1 

K the nondimensionalized variables x' = — , y' = -f-, a' = — , u* = — , t = - , 

d d a a a a d 

where d is the diameter of the jet and a a is the ambient speed of sound, are used, 

then equations (6), (7), (13), and (14) become the nondimensional set 


dx’. , , 

dt’ U + l * 


( 20 ) 


d£_ = 

dt’ 


ma' 


dZ_ _ 
dt’ 


dm 

dt' 


9u' 9a’ ,2 9u' 

V 9x' 9x’ 


, m 9u’ . _ 9a' 7 

cm — — r + m 7 — t 
l 9x' 8x' 


9y* 

2 £u]_ 
9y* 



with the initial conditions 


( 21 ) 

( 22 ) 

(23) 


x'(0) = y'(0) = 0 Z(0) = cos 6 m(0) = sin 9 

In order to approach these equations numerically, it is necessary that the velocity 
and temperature profiles be known. In an experimental situation, these profiles could be 
determined with probes. However, in the examples in this analysis, particular mathe- 
matical expressions were used for these profiles. These expressions were chosen 
because they were convenient computationally and illustrative of conditions occurring in 
an actual jet. However, the qualitative results are in no way dependent upon the partic- 
ular expressions chosen. 


Example 1 


Velocity and temperature profiles .- For the velocity profile, it is assumed that 


u’(x',y') 


M a 

^1 + /3x’ 



(24) 


This expression was employed because it is similar to that found in a normal spreading 
jet. Here, M a is the Mach number of the flow in the center of the jet-exit plane with 
respect to the ambient speed of sound. The constant a determines the shape of the 


9 



profile and the constant ft determines the rate at which the jet spreads. This velocity 
profile is shown in figure 3 for several values of /3x'. Note that the radius of the jet 
pipe corresponds to y’ = 1, that is, ay’ = a. 

A measure of the rate at which the jet spreads may be obtained by noting that 
99 percent of the velocity variation will occur in -3a < y' < 3a, where 


a = 



(25) 


Thus, the rate of spread is related to da/dx’. The spread angle <p is defined by 


tan 



3/3 

2 / 2 a 


(26) 


For the temperature profile, two variations will be considered, one having a high- 
temperature core and the other a low-temperature core. Comparison of these will allow 
evaluation of the effects of altering the temperature profile on the refraction of sound by 
the jet. 


Variation 1.- High -temperature core: 


T H (x',y’) = T + - A exp 
11 a \/l + yx' 


1 + yx') 


(27) 


This profile is shaped like the velocity profile obtained from equation (24). Thus, the 
maximum temperature is on the jet axis, the state which normally occurs in a jet. 


Variation 2.- Low -temperature core 

2A 


T L (x\y f ) = T a + 


\fl + yx' 


exp 


! (y ' 2 + y* 2 ) 

1 + yx' 


26 2 y’y* 
1 + yx' 


cosh 


(28) 


This profile is shown in figure 4 for several values on yx\ At the exit plane, it reaches 
a maximum at some point removed from the jet axis. 

For the purpose of comparison, the space average temperature deviation at corre- 
sponding cross sections has been required to be equal for both low- and high -temper ature - 
core profiles. This condition 



yields 


5 2 = 26 1 

Ray solutions .- In figures 5 and 6, the ray solutions y'(x') of equations (20) 
to (23) are shown for 0 fll 90° in steps of 6° for the parametric set 

M a = 0.887 qi 2 =1.6 /3=0.0 A = 2.76T a 

y = 0.0 6j 2 = 3.0 6 2 2 = 12.0 y* = 0.8 


This parametric set corresponds to a jet with a maximum exit velocity of 1000 ft/sec 
(305 m/sec) and a maximum temperature of 1500° F (1089° K). Because and y 
are both zero, this jet does not spread. The velocity and low -temperature -core profiles 
for this jet are shown in figures 3(a) and 4(a), respectively. In figure 5, where the high- 
temperature -core profile was employed, only the singular ray for 0=0° was channeled. 
However, in figure 6, the low-temperature -core profile channeled all rays for which 
0 ^ 42°. The channeled rays are periodic and never escape from the channel. 

Since this jet does not spread, it must follow the analysis developed previously for 
an axially uniform jet. The existence of a y c such that 




±1 


(29) 


was determined to be a necessary and sufficient condition for channeling. In figure 7, the 
points y' for which this condition is met are shown as functions of the initial angle 9 

for the parametric set previously mentioned. Thus, the jet will actually channel all rays 
for which 9 S 43°. 

Finally, in figures 8 and 9, the same parametric set has been employed except that 
/3 = 0.188 and y = 0.2. The use of these values causes the jet to spread in the axial 
direction. Here the spread angle <p, given by equation (26), is approximately 10°. The 
a, b, c, d, and e parts of figures 3 and 4 show relevant profiles which are approximately 
0.0, 2.5, 5.0, 7.5, and 10.0 jet radii from the jet -exit plane, respectively. In figure 8, 
where the high -temperature -core profile was utilized, again only the 0=0° ray was 
channeled. In figure 9, it can be seen that the low -temperature -core profile again initially 
channeled all rays for which 0 S 42°. Thus, it appears that equation (29) is approxi- 
mately valid for spreading jets as well. However, the rays for the spreading jet are not 
periodic and ultimately escape from the channel. 
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Example 2 


Velocity and sound-speed profiles.- In this example, the velocity and speed of sound 
are taken as polynomials since this is a simple method of describing experimental data. 
Knowledge of the sound speed is equivalent to knowledge of the temperature profile since 
they are directly related by equation (3). For the velocity profile, it is assumed that 


u’(x’,y’) = 


| M a + V' 2 + b 2 y ' 4 + b 3 y ' 6 


(-1 S y' * 1)1 


(Otherwis e)J 


and for the speed of sound, 


a’(x',y’) = 


|1 + Cjy' 2 + c 2 y’ 4 + c 3 y’ 6 


(-1 £y’ si)| 
(Otherwise) 


(30) 


(31) 


This jet does not spread, but spreading could be allowed by supposing the b^’s, c^’s, 

and the bounds on y’ to be functions of x’ . 

Ray solutions .- In figure 10, the polynomial velocity and sound-speed profiles are 
shown for the parametric set 


M a = 0.89 

b x = -5.96 

b 2 = 12.14 

b 3 = -7.07 

c^ = 6.87 

c 2 = -13.32 

Cg = 6.45 



These parameters correspond to a jet with maximum velocity 1000 ft/sec (305 m/sec), 
which attains a maximum temperature of 1500° F (1089° K) in an off-axis position. 

In figures 11 to 13, the ray solutions produced by these profiles are shown for var- 
ious initial conditions. Figure 11 depicts the rays for -90° = 0 = 90° in steps of 6° for 
the initial condition x'(0) = y’(0) = 0. All rays for which -24° 2 6 £ 24° are channeled. 
In figure 12, the ray solutions are displayed for the initial condition x'(0) = 0, 
y'(0) = 0.316. For these conditions, only those rays for which -12° § Q < 12° were 
channeled. Figure 13 shows the ray solutions when x'(0) = 0, y’(0) = 0.706. No rays 
are channeled under these conditions. Finally, in figure 14, the singular values y’ are 
shown as a function of 6 for the initial conditions in which rays were channeled. For 
y'(0) = 0, all rays for which -29° = 9 = 29° would be channeled and for y'(0) = 0.316, 
all rays for which -16° = 9 = 16° would be channeled. 
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CONCLUDING REMARKS 

In this analysis, the techniques of geometrical acoustics have been employed to 
study the effect of variation of the jet-temperature profile upon propagation of sound in 
the jet. It has been shown that the temperature profile of the jet can be shaped to cause 
the jet flow to act as a channel for the sound. Further, a simple analytical expression 
which predicts the occurrence of such phenomena in a nonspreading jet has been obtained. 
This expression also seems to be approximately valid for a spreading jet. 

From this analysis, the conclusion that jet noise may be favorably affected by 
varying the jet -temperature profile can be drawn. However, it should be mentioned that 
no allowance for changes in the production of noise by the jet has been included; that is, 
the mixing process has been considered to be independent of the temperature gradient. 

It is conceivable that altering the jet -temperature profile could increase the noise of the 
jet and ultimately negate any favorable effects of channeling the sound. Further, the 
important question of the thrust efficiency of such an engine has not been evaluated. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., January 30, 1970. 
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